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INTRODUCTION 

The potential for the use of metal matrix composite materials in structural 
applications is, in many respects, proportional to such factors as the extent of 
material characterization, the practical experience with the material, and the 
applicability of computational methods. The content of this paper is directed 
toward adding to the general knowledge of these factors with particular emphasis 
on the prediction of thermal stresses in a hybrid metal matrix composite spar 
cap element (references 1 and 2). The need to critically evaluate computational 
procedures has increased in recent years because the cost of laboratory testing 
and test specimens has resulted in a deficiency of test data providing a basic 
comparison of calculated and measured thermal stresses (reference 3 ). The wide 
spectrum of configuration/material options has also contributed to this 
deficiency. 

A hybrid metal matrix composite upper surface wing spar cap element 
(reference 2) composed of titanium which was selectively reinforced with boron- 
carbide-coated boron filaments in an aluminum matrix was used as a specimen to 



compare computed thermal stresses with measured values obtained from a thermal 
cycling experiment. The instrumented test component was heated in an oven and 
the thermal strains resulting from dissimilar material effects were measured. 

The metal matrix composite portion of the spar is an orthotropic unidirectionally 
reinforced fibrous composite material. A finite element NASTRAN (reference 4) 
model with appropriate anisotropic elements was developed to provide 
computational data for comparative purposes. Rules of mixtures were utilized 
to estimate anisotropic elastic and expansive properties and these results were 
incorporated into the structural model for comparison with experimental results. 
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SYMBOLS 
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Young f s modulus, Pa (lbs/in ) 
non-linear mathematical representations 
shear modulus, Pa (lbs/in ) 
general linear constant 
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stiffness constant. Pa (lbs/in ) 
temperature, incremental temperature, K (F° ) 
volume fraction 

coefficient of linear expansion, m/mK (in/inF c ) 
strain, m/m ( in/in ) 
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SUBSCRIPTS 
refers to apparent 
refers to corrected value 
refers to fiber 
arbitrary numbers 
refers to matrix 
refers to standard 
refers to metal matrix composite 
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DESCRIPTION OF TEST 


The origin of the hybrid titanium/metal matrix composite test specimen 
and the associated reasoning behind the configuration is well documented in 
reference 2. The specimen is configured to be associated with the upper spar 
cap area as is depicted in figure 1. The fibers are oriented unidirectionally 
along the longitudinal axis of the test specimen. The 1.09 meter (36.0 inch) 
specimen is shown in the photograph of figure 2. This figure also shows end 
attachments which were installed to accommodate tensile load tests not discussed 
in this paper. A dimensional sketch of a segment near the middle of the 
specimen which locates the principal strain gages used for comparison with the 
structural model is shown in figure 3. Thermocouples and strain gages are 
located as shown in figure 4. The presence of additional strain gages which 
were utilized in axial load tests not presented in this paper should be noted. 
Foil-type resistance strain gages and Chromel-alumel thermocouples were 
utilized for this test. 

The data presented in this paper was obtained from a single thermal cycle 
of a specimen that was undergoing a 200 cycle thermal exposure test in an oven. 
The nature of the thermal cycle is depicted in the time-history of temperature 
presented in figure 5. A cycle of thermal strain resulting from the elevated 
temperature of the specimen is also shown in figure 5. The basic test plan 
for each individual cycle was to raise the temperature from ambient to a value 
of 47 8 K (400 F" ), maintain the temperature at this value for thirty minutes, 
and then return the specimen temperature to ambient within a total cycle time 
of two hours. A photograph of the cycling oven is shown in figure 6. 

Because of the anisotropic character of the metal matrix composite, a 
major problem associated with this laboratory experiment was the appropriate 
selection of strain gages. The temperature compensation/apparent strain 
relations were critical since the thermal expansion properties in the three 
major material axes were not known precisely (this was not a problem with 
the titanium part of the specimen since the Ti-6A1-4V alloy has thermal 
expansion properties that are well documented). Strain gages mounted on the 
metal matrix composite transverse to the fibers (2-direction, see figure 3) 
were compensated for the aluminum alloy of the matrix. The strain gages on 
the metal matrix composite mounted parallel to the fibers ( 1-direction, see 
figure 3) were compensated for molybdenum since this readily available strain 
gage had a value which most closely represented what was thought to be the 
thermal expansion properties in the fiber direction. 

The precise acquisition of laboratory strain gage test data requires more 
than presumptions about the thermal expansion properties. Therefore, two 
strain gages were mounted on a small coupon of the metal matrix composite 
parallel to the fibers and two more were mounted transverse to the fibers. 
Precise apparent strain curves were obtained by heating the coupons very slowly 
in an oven and recording the resulting strain. This data was used to correct 
the test data obtained from the hybrid titanium/metal matrix composite specimen 
for apparent strain. This data was also used to establish values for the 
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thermal coefficient of expansion which will be discussed in detail in the 
following section. 


COEFFICIENT OF THERMAL EXPANSION 

The thermal expansion coefficients for the metal matrix composite were 
determined from the strain gaged coupon seen in figure 7. The precisely known 
thermal expansion properties of the strain gages were used as the standard. The 
specimen with attached strain gages was heated slowly in an oven up to the test 
temperature of 478 K (400 F s ). One pair of strain gages ( 1-direction) was 
compensated thermally for molybdenum and the other pair (2-direction) was comp- 
ensated for aluminum. Another pair of strain gages compensated for molybdenum 
were oriented in the 2-direction on another specimen not shown. This was done 
because the values in the 2-direction were quite uncertain in the mind of the 
author. The thermal expansion properties were deduced from deviations of the 
strain gages on the metal matrix composite from the strain gage manufacturers 
predicted apparent strain curves. 

Deducing the coefficient of thermal expansion of a material from attached 
strain gages is relatively straightforward since an accurate standard is 
inherent (the compensation value of the strain gage). An apparent strain curve 
(reference 5) is required and provided with strain gages because the coefficient 
of thermal expansion is not truly linear over the usable temperature range. 

If a strain gage of known compensation and known apparent strain is mounted 
to a material of unknown thermal expansion properties, then the thermal 
expansion properties of that material can be deduced from the deviation from 
the apparent strain curve using an oven test. The basic equation defining 
thermal displacement may be expressed as: 

otj^AT = a s AT + f(T) - g(T) (l) 

where the first term on the right hand side represents the linear portion 
of the standard material, the second term represents the non-linear deviation 
from the linear presumption of the first term, and the third term represents 
the deviation from the first two terms because the material on which the strain 
gage is mounted is different from the standard. If the quantity 

f ( T ) - g(T) (2) 

is determined by experiment, then a linear approximation can be made and a 
correction term may be identified as 

a AT - f ( T ) - g(T) 

hence 

V T = (a s + a c )AT 

or 


(3) 

(4) 
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(5) 


M 


a 
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The linear coefficient of thermal expansion in the fiber direction ( 1-axis) 
of the metal matrix composite was determined from the standard and corrective 


values: 

“M = “s + “c = • 0000 °54 + .0000009 = .0000063 

111 (.0000030) (.0000005) (.0000035) (6) 


The expansion coefficient in the transverse direction (2-direction) was 
determined using two strain gages of largely differing coefficients of thermal 
expansion. Hence, there were two sets of data used to determine the values. 
The test results indicate slightly differing values: 



a + a = .0000234 
2 2 (. 0000130 ) 


.0000036 
( . 0000020 ) 


.0000198 
( . 0000110 ) 


(7) 


and 



= .0000054 + .0000135 

(.0000030) (.0000075) 


.0000189 
( .0000105) 


( 8 ) 


A nominal value of .00001935 m/mK ( .00001075 in/inF" ) was used for the 
thermal expansion coefficient in the transverse (2-axis) direction of the metal 
matrix composite. 

A sample of the data used in determining the coefficients of thermal 
expansion is shown in figure 8. The manufacturers apparent strain curve, i.e., 
the f(T) term in equation (l), is represented by the solid lines. The deviation 
from the manufacturers apparent strain curve is represented by the difference 
between the line and the circular symbols. This deviation constitutes the g(T) 
term in equation (l). It should be noted that the deviation is cumulative strain 
and must be divided by the temperature change to arrive at a linear approximation 
of the change in coefficient of thermal expansion. 


STRESS-STRAIN RELATIONS 


The stress-strain relations in terms of stress components as a function 
of material properties and strain components is: 


Q. .e . 
ij J 


(i,j = 1, 2, 6) 


(9) 


Hence, for an orthotropic, unidirectional reinforced fibrous composite 
material : 
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( 10 ) 



The stiffness constants are as follows: 
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The reciprocity relation discussed in reference 6 requires Q..„ = Q 21 


because : 


21 


12 


( 12 ) 
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The basic elastic properties are needed as input for the NASTRAN model 
since the MAT2 (reference 4) inputs require the diagonal half of the material 
property matrix of equation 9: 


Q 


11 


'12 

Q 16 

'22 

Q 26 


Q 66- 


(13) 


The relationships of equations 10 and 11 are also needed to reduce the 
measured strain data acquired on the metal matrix composite portion of the 
structure. 


LAMINA ELASTIC PROPERTIES 


The accurate knowledge of elastic properties in composite materials is 
critical to their safe and effective use. Frequently, insufficient or incomplete 
information is available to develop full confidence in material characterizations 
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and assumptions. The hybrid metal matrix composite spar which is the subject 
of this paper is no exception with respect to material properties. The basic 
elastic properties of the hybrid metal matrix composite spar constituents were 
obtained from the combination of known constituent material properties and 
rules of mixtures theory. 


Young T s modulus in the direction parallel to the continuous fibers 
( 1-direction ) can be presumed from the rule of mixtures (references 6 and 7): 


E n = E„V„ + E V 
1 ii mm 


(14) 


The volume ratio and V m were determined to be 0.5 from photomicrographs 
(reference 2). Based on constituent material properties, the value of was 
determined to be 242.0 x 10^ Pa (35.1 x 10 G lbs/in^). 


The elastic modulus transverse to the fiber direction (2-direction) was 
computed from the formula: 


E„E 
f m 


Vf + V f E m 


(15) 


Q 2 

The value of E^ was determined to be 110.9 x 10 Pa (16.1 x 10 lbs/in ). 
The shear modulus for the metal matrix composite, G , was estimated from 


the equation ( reference 6 ) : 


12 


w 1 + v » v 1 - h » 

Gil - V.) * Gil + T.) 


The relationships, 


(16) 


G V 


2(1 + v f ) 


G ~ 
m 


m 


2(1 + v ) 
m 


(17) 


yield the individual values for the fiber and the matrix respectively. The 
magnitude of G ^ was estimated to be 59.61 x 10 G Pa (8.65 x 10^ lbs/in^). 

The rules of mixtures for Poisson's ratio and the individual fiber and 
matrix properties can be used to arrive at a value for 


v 


12 


V v 
m m 


Vf 


( 18 ) 
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The values V = V„ = 0.5, v = .33, and v. = .13 lead to the value for v 
m f m 1 

which is computed to be = . 23 . 

Since the compliance matrix is symmetric, then: 


12 


v 12 _ V 21 


v 12 E 2 


21 


- .11 


(19) 


Summarizing the basic elastic properties used for analysis: 

E x = 241.9 x 10 9 Pa (35.1 x 10 6 lbs/in 2 ) 

E 2 = 110.9 x 10 9 Pa (16.1 x 10 6 lbs/in 2 ) 

G 12 = 59,6 x 1C)9 Pa ^ 8#65 x 1C)6 lbs / in2 ) 

V 12 = * 23 
V 21 " ,:L1 

The matrix of stiffness constants developed in equations 10 and 11 can be 
computed from the previous basic elastic properties: 


or in English units: 


178. 6 
19.4 
0 

25.90 

2.81 

0 


19.4 

84.1 

0 

2.81 

12.20 

0 


0 

0 

59.6 

0 

0 

8.65 


x 10' Pa 


x 10 G lbs/in 2 


The use of the matrix of stiffness constants is discussed in the next 
section where the structural model used for predictions is presented. 


STRUCTURAL MODEL 

A finite element model of the hybrid metal matrix composite spar element 
was developed using NASTRAN (reference 4). A sixty plate element (CQUAD2), 
sixty-eight grid point model was developed for a symmetrical half of the spar. 
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The schematic of the model is shown in figure 9. The model was specifically 
contructed to correlate with strain gages on the actual specimen. The metal 
matrix composite material is represented by twelve plate elements in the 1-2 
plane. All elements in the 1-3 plane are titanium. The only loading input to 
the model was a temperature change at the grid points consisting of the 
difference between room temperature and the maximum oven temperature of 478K 
(400 F c ). The stiffness constants developed in the previous section are used 
as input to NASTRAN MAT2 cards in the bulk data deck. 


RESULTS AND DISCUSSION 

The basic thrust of this paper has been directed toward evaluating 
analytical methods to predict thermal stresses in a reinforced fibrous 
composite material. Since the test structure was a relatively simple component, 
the structural modeling turned out to be relatively straightforward while the 
definition of the material properties of the metal matrix composite was 
deceptively complicated. The degree of correlation of predicted thermal 
stresses with measured values will tend to reflect how well the metal matrix 
composite material was characterized rather than how well the overall spar 
was modeled. 

Thermal strain measurements were taken on the hybrid metal matrix 
composite spar at the time indicated in figure 5 at locations corresponding 
to elements A, B, and C (figures 9 and 10). The two elements in the 1-3 plane, 
elements A and C, represent titanium parts of the spar while the element in 
the 1-2 plane, element B, represents the metal matrix composite part of the 
spar. The circular symbols shown in figure 10 represent thermal stresses 
obtained by reducing strain data obtained from strain gages mounted on the 
spar during the thermal cycling. Values for the elastic properties used to 
reduce the strain gage data are presented in the test measurement row in 
Table I. The solid lines shown in figure 10 represent computed thermal stresses 
based on the structural model shown in figure 9. The anisotropic elastic 
properties are determined in the Lamina Elastic Properties Section and are 
presented in Table I. 

The correlation between predicted and measured values of thermal stresses 
should be viewed with the recognition that the metal matrix composite and its 
constituents have not been acceptably characterized. This has great impact 
on the results since Young ! s modulus and Poisson’s ratio are used in the 
reduction of the strain data and in the computations developed with the 
structural model. An anisotropic structure is particularly sensitive to 
variations in material properties. These effects have apparently led to some 
large errors in correlating measured and calculated thermal stresses. 

Data of a different computational nature is presented also in Table I 
that provides insight into the role that the anisotrophy and the material 
property inaccuracy of the metal matrix composite play in the thermal stress 
computation. The basic structural model was rerun for the same temperature 
case, however, the material properties for the metal matrix composite were 
presumed to be isotropic. In one case the metal matrix composite elements 
were presumed to be isotropic with the elastic and expansional properties 
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corresponding to those in the 1-direction. The other case was run using only 
the 2-direction properties (refer back to the Lamina Elastic Properties Section 
and the Coefficient of Thermal Expansion Section). The resulting correlation 
is seen (Table I) to degrade significantly from the anisotropic to the isotropic 
analysis. Since the change that has occurred is a material property change, it 
must be concluded that failure to account for anisotrophy and failure to 
accurately characterize the material are most likely dominant factors in the 
correlation of measured and computed thermal stresses of the hybrid metal matrix 
composite spar presented in this paper. This illustrates clearly the importance 
and necessity of a well defined, meaningful materials characterization prior 
to concept or computational tool assessments. 


CONCLUDING REMARKS 

A hybrid spar of titanium with an integrally brazed composite consisting 
of an aluminum matrix reinforced with boron-carbide-coated boron fibers was 
heated in an oven and the resulting thermal stresses were measured. Uniform 
heating of the spar in an oven resulted in thermal stresses arising from the 
effects of dissimilar materials and the anisotrophy of the metal matrix 
composite. The thermal stresses were calculated from a finite element 
structural model using anisotropic material properties deduced from constituent 
properties and rules of mixtures. 

Comparisons of calculated thermal stresses with measured thermal stresses 
on the spar are presented. It was shown that failure to account for 
anisotrophy in the metal matrix composite elements would result in large 
errors in correlating measured and calculated thermal stresses. It was also 
concluded that very strong material characterization efforts are required to 
predict thermal stresses well in anisotropic composite structures. 


NASA Awes Research Center 
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February 28 y 1985 


REFERENCES 

1. Sakata, I.F.j and Davis, G.W. : Evaluation of Structural Design Concepts 
for an Arrow-Wing Supersonic Cruise Aircraft. NASA CR-2667, April, 1977. 

2. Sakata, I.F.j Davis, G.W.; and Cox, J.M. : Advanced Hybrid Structural 
Concepts Development - Technology Study for Advanced Supersonic Cruise 
Vehicles. NASA CR-3728, 1983. 

3. Wagner, R.E.j Stultz, J.W.j Lee, R.C.; and Wilkinson, R.D. : Development 
of Design Criteria and Procedures to Account for the Effects of Elevated 
Temperatures Upon Missile Airframes - Tasks I, II, and III. AFWAL-TR-82- 
3099, December, 1982. 


10 



4. NASA-COSMIC: NASTRAN Users Manual. NASA SP-22l(05), 1978. 

5. Jenkins, Jerald M. : A Study of the Effect of Apparent Strain on Thermal 
Stress Measurement for Two Types of Elevated Temperature Strain Gages. 

NASA TM-84904, February, 1983. 

6. Whitney, James W. ; Daniel, Isaac,; and Pipes, R. Byron: Experimental 
Mechanics of Fiber Reinforced Composite Materials; Society for Experimental 
Stress Analysis, Monograph No. 4, 1982. 

7. Jones, Robert M. : Mechanics of Composite Materials. McGraw-Hill Book 
Company, 1975. 

8. Jenkins, Jerald M. ; Schuster, Lawrence S. ; and Carter, Alan L. : Correlation 
of Predicted and Measured Thermal Stresses on a Truss-Type Aircraft 
Structure. NASA TM-72857, 1978. 

9. Jenkins, Jerald M. : Correlation of Predicted and Measured Thermal Stresses 
on an Advanced Aircraft Structure with Similar Materials. NASA TM-72862, 
1979. 

10. Jenkins, Jerald M. : Correlation of Predicted and Measured Thermal Stresses 
on an Advanced Aircraft Structure with Dissimilar Materials, NASA TM-72865, 
1979. 


11 



TABLE I. - SUMMARY OF MATERIAL PROPERTIES AND THERMAL STRESSES 


(a) Material properties 


ELASTIC 

PRESUMPTION 

TITANIUM 

METAL MATRIX COMPOSITE 

a 

E 

V 

“1 

°2 

E 1 

E 2 

v 12 

V 21 

Isotropic 

9.5 x 10-6 

100.0 X 10 9 

• 31 

19.4 x 10-6 

19.4 x 10-6 

82.0 x 10 9 

82.0 x 10 9 

.11 

.11 

(1-direction 

(5.3 x 10-6) 

(14.5 x IQ®) 

(.31) 

(10.8 x 10-6) 

(10.8 x 10-6) 

(11.9 x 106 ) 

(11.9 x 10®) 

(.id 

(.11) 

properties) 










Isotropic 

9.5 x 10-6 

cn 

o 

X 

o 

o 

o 

• 31 

6.3 x 10-6 

6.3 x 10-6 

173.7 x 109 

173.7 x 10 9 

.23 

.23 

( 2-di rection 

(5.3 x 10-6) 

(14.5 x 106 ) 

(.31) 

(3.5 x 10-6) 

(3.5 x 10-6) 

(25.2 x 106) 

(25.2 x 106 ) 

(.23) 

(.23) 

properties) 










Anisotropic 

9.5 x 10-6 

o 

o 

o 

X 

o 

.31 

6.3 x 10*6 

19.4 x 10-6 

241.9 x 10 9 

110.9 x 10 9 

.23 

.11 


(5.3 x 10-6) 

(14.5 x 106 ) 

(.31) 

(3.5 x 10-6) 

(10.8 x 10-6) 

( 35.1 x 106 ) 

(16.1 x 10 6 ) 

( • 23) 

(.11) 

Test Measurement 

NA 

100.0 x 10 9 

.31 

NA 

NA 

241.9 x 10 9 

110.9 x 10 9 

.23 

.11 


NA 

(14.5 x 10 6 ) 

(.31) 

NA 

NA 

(35.1 x 106 ) 

(16.1 x 106 ) 

(.23) 

(.11) 

SI Units 

m/mK 

Pa 


m/mK 

m/mK 

Pa 

Pa 

none 

none 

(English Units) 

(in/inF°) 

(lbs/in 2 ) 


(in/inF°) 

(in/inF°) 

(1 bs/i n 2 ) 

(lbs/in 2 ) 

none 

none 


(b) Thermal stress 


ELASTIC 

ELEMENT A 

ELEMENT B 

ELEMENT C 

PRESUMPTION 

°i 

°i 

°2 

°i 

Isotropic 
( 1-di rection 
properties) 

49.5 x 10 6 
(7185) 

76.5 x 10 6 
(-11100) 

40.0 x 10 6 
(-5810) 

62.5 x 10 6 
(9060) 

Isotropic 

(2-direction 

properties) 

-27.2 x 10 6 
(-3940) 

43.2 x 106 
(6260) 

24.8 x 106 
(3600) 

-36.4 x 10 6 
(-5280) 

Anisotropic 

-28.8 x 10 6 

(-4174) 

26.8 x 10 6 

(3880) 

45.2 x 10 6 
(-6560) 

-39.5 x 10 6 
(-5 7 30) 

Test Measurement 

-24.5 x 106 
(-3550) 

45.2 x 106 
(6550) 

-70.9 x 106 
(-10280) 

-26.9 x 10 6 
(-3900) 

SI Units 
(English Units) 

Pa 

( 1 bs/i ) 

Pa 

(1 bs/i ) 

Pa 

(lbs/in?) 

Pa 

(lbs/in?) 
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Figure 7 » Metal matrix composite coupon used to determine 
thermal expansion properties . 
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Figure 8 * Strain gage apparent strain deviation from manu- 
facturer prediction for metal matrix composite coupons • 






Figure 9 . Finite element model of hybrid 
metal matrix composite spar cap . 



Figure 10. Comparison of measured and calculated ther- 
mal stresses • 
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